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Abstract 

We develop a systematic way for constructing bispectral algebras of commut- 
ing ordinary differential operators of any rank N . It combines and unifies the 
ideas of Duistermaat-Griinbaum and Wilson. Our construction is completely 
algorithmic and enables us to obtain all previously known classes or individual 
examples of bispectral operators. The method also provides new broad families 
of bispectral algebras which may help to penetrate deeper into the problem. 
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Introduction 



In this paper we reconsider the bispectral problem. As stated in [DG], it asks for 
which ordinary differential operators L(x, d x ) there exists a family of eigenfunctions 
*$>(x,z) that are also eigenfunctions for another differential operator A(z, d z ) but this 
time in the "spectral parameter" z, to wit 



L(x,d x )V(x,z) = f(z)¥(x,z), 
A(z,d z )^(x,z) = Q(x)^(x,z) 



(0.1) 
(0.2) 



for some functions f(z),@(x). Both operators L and A are called bispectral. 

This problem first appeared in Gl ] in connection with "limited angle tomogra- 
phy" (see also |G2, |G3| , DG]). Later it turned out to be related with several, seem- 
ingly far from it, topics and in particular, with soliton mathematics. To be more 
specific, we have to mention the deep connection with some very actively developing 
areas of research in mathematics and theoretical physics like Calogero-Moser parti- 
cle system | W2, Kf (see also Q), additional symmetries of KdV and KP hierarchies 



BHY4], representation theory of Wi +00 -algebra [BHY4], etc. These studies 



not only revealed the rich mathematical structure behind the bispectral problem, 



but also (if we use a remark by G. Wilson | W2 |) "deepened the mystery" around it. 
Thus, not only applications, but also purely mathematical questions motivated the 
great activity in the past few years in the bispectral problem. 

In the present paper we construct new families of bispectral operators. In order 
to explain better our contribution, we need to review some of the achievements in 
the subject. 

The first general result in the direction of classifying bispectral operators belongs 
to J. J. Duistermaat and F. A. Griinbaum |DG|| . They determined all second order 
operators L admitting an operator A such that the pair (L, A) solves the bispectral 



problem (0.1, D.2). Their answer is as follows. If we write the operator L in the 
standard Schrodinger form 

d2 

+ u(x), 



L 



dx 2 



the bispectral potentials u(x) are given (up to translations and rescalings of x and 
z) by the following list: 



u{x) 
u(x) 



cx 



c e C 



(Airy); 
(Bessel); 



(0.3) 
(0.4) 



u(x), which can be obtained by finitely many rational Darboux 

transformations from u(x) = 0; (0-5) 

u(x), which can be obtained by finitely many rational Darboux 

1 

Ax 2 ' 



transformations from u(x) 



(0.6) 



The family (|0.5| ) has previously appeared in [ AMM , AM] and is known as "ratio- 
nal solutions of KdV" . They can be obtained also by applying "higher KdV flows" 
to potentials Vk{x) = k(k + l)x~ 2 , k £ N. 
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The second family ( p.6| ) was interpreted by F. Magri and J. Zubelli |MZ| | as 
potentials invariant under the flows of the "master symmetries" or Virasoro flows. 
Besides the classification of the bispectral operators by their order, another 



scheme has been suggested in [DG] and used in [Wl]. Below we explain it in a 
general context as it will be used throughout this paper. One may consider an 
operator L(x,d x ) as an element of a maximal algebra A of commuting ordinary dif- 



ferential operators BC]. Following G. Wilson [Wl], we call such an algebra bispectral 
if there exists a joint eigenfunction ^f(x, z) for the operators L in A that satisfies 
also equation (|0.2| ). The dimension of the space of eigenfunctions ^f(x,z) is called 



rank of the commutative algebra A (see e.g. [ KrN|] ). This number coincides with 



the greatest common divisor of the orders of the operators in A. For example, the 
operators with potentials (p.5|) belong to rank 1 algebras and those with potentials 



(fjg , Qj, gj) to rank 2 algebras 

All rank 1 maximal bispectral algebras were recently found by G. Wilson [Wl] 
These algebras do not necessarily contain an operator of order two. 



The methods of the above mentioned papers [DG| and [Wl| may seem quite 



different. Indeed, while in [ DG | the "rational" Darboux transformations play a 
decisive role, G. Wilson Wl ] uses planes in Sato's Grassmannian obtained from the 
standard H + = span{z fc }fc>g by imposing a number of conditions on it. One of our 
main observations is that both methods, appropriately modified, can be looked upon 
as the two sides of one general theory. 

From this new point of view in the present paper we construct nontrivial maximal 
bispectral algebras of any rank N, thus extending the results from [DG, Wl]. For 



example, for any positive integer k we obtain bispectral algebras of rank N with the 
lowest order of the operators equal to kN. Our method allows us to obtain all classes 
and single examples of bispectral operators known to us by a unique method. At the 
same time we suggest an effective procedure for constructing bispectral operators, 
despite the fact that the theory involves highly transcendental functions like Airy or 
Bessel ones. The point is that the latter are used in the proofs while the algorithm 
given at the end of Sect. 3 performs arithmetic operations and differentiations only 
on explicit rational functions. 

In the rest of the introduction we describe in more detail the main results of the 
paper together with some of the ideas behind them. 

The framework of our construction is Sato's theory of KP-hierarchy |S], [DJKMl 
|W| , IvMfl . In particular, our eigenfunctions are Baker or wave functions ^>y(x,z) 



corresponding to planes V in Sato's Grassmannian Gr and our algebras of commut- 
ing differential operators are the spectral algebras Ay. We obtain our bispectral 
algebras by applying a version of Darboux transformations, introduced in our pre- 
vious paper ( |BHY3U , on specific wave functions which we call Bessel (and Airy) wave 
functions (see Sect. 1 and 4). As both notions are fundamental for the present paper 
we hold the attention of the reader on them. Bessel wave functions are the simplest 
functions which solve the bispectral problem (see Jz| where they were introduced 
and 0). They can be defined as follows. For (3 = (J3i,..., (3 N ) G ^p{x, z) is the 
unique wave function satisfying 

xd x ^p(x,z) = zd z ^p(x,z) 
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(i.e. typ(x, z) depends only on xz) and 

Lp(x,d x )^fs(x,z) =z N ^f3(x,z), 

where Lp(x, d x ) = x~ N (xd x — Pi) ■ ■ ■ (xd x — (3n) is the Bessel operator. Obviously, 
the above equations lead to 

Lp{z,d z )^p{x,z) =x N ^ (3 (x,z). 

Similarly, for a = (ao,«i, • • • , ajv-l) £ C^ -1 consider the (generalized) Airy wave 
function (see ]KS| , |Dij| ] ) satisfying: 

(d? + J2 a ^' 1 - Z ) = zN ^a{x, Z). 

It depends only on oqx + z N and again gives a simple solution to the bispectral 
problem. The Airy case is in many respects similar to the Bessel one. As we find 
the latter case richer in properties, we pay more attention to it, contenting ourselves 
only with a sketch of the former. 



Classically, a Darboux transformation [BC, Da| of a differential operator L, pre- 
sented as a product L = QP, is defined by exchanging the places of the factors, 
i.e. L = PQ. Obviously, if 9(x, A) is an eigenfunction of L, i.e. L(x, d x )^(x, A) = 
\ty(x, A), then P^/(x, A) is an eigenfunction of L. 

Here we introduce Darboux transformations not only on individual operators 
but also on the entire spectral algebra corresponding to a Bessel (or Airy) plane. In 
other words, we apply them on operators L which are polynomials h(Lp) of Bessel (or 
Airy) operators. These transformations may be considered as Backlund-Darboux 



transformations on the corresponding wave functions [ AvM | . Such Darboux trans- 
formation is completely determined by a choice of a Z^-invariant operator P(x,d x ) 
with rational coefficients normalized appropriately by a factor g~ 1 (z) to ensure that 

* w (x,z) = ——P{x,d x )^p{x,z) 
9{z) 

is a wave function. We call ~$fyy(x, z) (respectively W) a polynomial Darboux trans- 
formation of ^ p (respectively Vn). The definition of polynomial Darboux trans- 
formations of Airy planes is similar to that in the Bessel case with only minor 
modifications: P is not necessarily Z^y-invariant and g(z) has to belong to Cfz^]. 
Thus we come to our main result. 

Theorem 0.1 If the wave function *&\y(x, z) is a polynomial Darboux transforma- 
tion of a Bessel or Airy wave function ^p{x,z), then it is a solution to the bispec- 
tral problem, i.e. there exist differential operators L(x,d x ), A(z,d z ) and functions 
f(z), e(x) such that |0]) and (Q hold. 

Note the difference between the classical definition and the definition introduced 



here. In contrast to [ DG | where the authors make a finite number of "rational" 
Darboux transformations, we perform only one polynomial Darboux transformation 
to achieve the same result. 
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Our definition of polynomial Darboux transformation is constructive as P(x, d x ) 
is determined by the finite dimensional space KerP. For this reason one can explicitly 
present at least one operator L € Aw) it can be given by Ph(Lp)P~ l . Usually it is 
of high order. But as it is only one element of the whole bispectral algebra Aw there 
can be eventually operators of a lower order. For example, the bispectral operators 
of [DG] are of order two. There is a simple procedure (see [ BHY3H ) to produce the 
entire bispectral algebra Aw of commuting differential operators. In addition, one 
can show that the spectral curve Spea4vK ( see e -g- [AMcD] for definition) is rational, 
unicursal and Z n -invariant. 

In the course of our work we have widely used important ideas introduced by 



G. Wilson [Wl]. Among them we mention first the idea of explicitly writing con- 
ditions on vectors of a plane V € Gr which define the new plane obtained by a 
Darboux transformation. Second is the notion of involutions on the Sato's Grass- 
mannian. In particular, we extend the bispectral involution b introduced in [Wl] to 
the manifolds of polynomial Darboux transformations. More precisely, we prove the 



following theorem, from which Theorem 0.1 is an obvious consequence. 



Theorem 0.2 (i) The bispectral involution is defined for planes W which are poly- 
nomial Darboux transformations of Bessel or Airy planes. 

(ii) The image bW of such a plane W is again a polynomial Darboux transfor- 
mation of the corresponding Bessel (respectively) Airy plane. 



Our main concern in the present paper is to prove Theorem 3.2. Our second 



goal is to provide explicit formulae and examples (see Sect. 5), which are not only 
an illustration of our method but also show the existence of new families of bispec- 
tral operators with particular properties. Some of them generalize directly the well 
known ones like the Duistermaat-Grunbaum's "even case" ( p.6| ) [DG|. Other fami- 
lies exhibit quite different properties from the well known examples. In this respect 
Sect. 5 has also the role to supply diverse experimental material for new insights 
into the theory of bispectral algebras. We draw the attention of the reader also 
to the explicit formulae for the action of the bispectral involution on an important 
class of Darboux transformations (which we call monomial) of Bessel operators. As 
a particular case, we obtain such formulae for all second order bispectral operators 
found in JDGJ. 

The class of monomial Darboux transformations has also other remarkable prop- 
erties, e.g. they are connected to representation theory of VFi+oo-algebra. We do 
not touch this matter here for lack of space. The interested reader can learn about 
it in |BHY4| . 

A natural question is if the operators found in this paper form the entire class 
of bispectral operators. The answer is negative as recently shown in [BHY5]. 

At the end for the reader's convenience we give a brief description of the or- 
ganization of the paper. Sect. 1 is intended only for reference. It reviews results 
connected with Sato's theory, which we need for the treatment of the bispectral prob- 
lem. Besides the general notions (see e.g. [^|, DJKM, |SW| ) we recall the involutions, 



introduced by G. Wilson Wl] and in particular, the bispectral involution. In Sect. 2 



we introduce our manifolds Gr^ of polynomial Darboux transformations of Bessel 
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planes. We give two equivalent definitions (Definition |2.5| and the one provided by 
the statement of Theorem |2.7|) . Sect. 3 contains our main results - Theorems |0.1| and 
[T^. for the Bessel case. Sect. 4 deals with the analogs of Sect. 2 and 3 for the Airy 
case (although in different order). The last Sect. 5 is devoted to explicit examples of 
bispectral operators, which have been studied in other papers | DG , [W2f| , as well as 
new families (which we have not seen elsewhere). The emphasis in Sect. 5 is rather 
on the simple algorithmic way of constructing bispectral operators (wave functions, 
etc.) than on the novelty of the examples. 

For readers who wish to see the main results as soon as possible we propose 
another plan of reading the paper. They can start with Sect. 2 and read it up to 



the statement of Theorem 2.7, returning to Sect. 1 for reference when needed. Then 



skipping the (technical) proof of Theorem 2.7, they can go Sect. 3. After that, taking 



for granted the proof of Theorem |3.2| , they can look at the examples of bispectral 
operators, originating from Bessel ones in Sect. 5. Thus they will have a complete 
picture of the results in the Bessel case, and having this experience, they can easily 
go through the Airy case. 

More detailed information about the material included in each section can be 
found in its beginning. 



The present paper is a part of our project on the bispectral problem BHY2]- 
TBHY5I1 . The main results contained here were announced at the conference of 



Geometry and Mathematical Physics, Zlatograd 95 (see [BHY1| 



After this paper was written, we got a paper [KR] where some of the results 
about the Airy case were obtained independently. 
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1 Preliminaries 

In this section we have collected results about Sato's theory, relevant to the bispectral 
problem. For reader's convenience we have divided the section into 4 subsections, 
whose titles, hopefully, give an idea of their content. The reader, who is acquainted 
with Sato's theory may even skip this section and return to it for references when 
needed. More detailed account of the material of the subsections can be found in 
their beginnings. 

1.1 Sato's Grassmannian and KP— hierarchy 

We shall recall some facts and notation from Sato's theory of KP-hierarchy needed 
in the paper. The survey below cannot be used as a systematic study. There are 
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several complete texts on Sato's theory, starting with the original papers of M. Sato 
and his collaborators DJKM] (see also [ 5W , 
Consider the space of formal series 



V = a,kV k Ofc = for k S> o| . 



Sato's Grassmannian Gr |S|, DJKM , [SW| ] consists of all subspaces (planes) W C V 
which have an admissible basis 

w k = Vk + ^WikVi, k = 0,1,2,... 

i<k 

In Sato's theory V is most often realized as the space of formal Laurent series in 
z^ 1 via v k = z k . The Baker (or wave function) ty\y(x, z) of the plane W contains the 
whole information about W as the vectors uik = d x ^w( x > z)\ x= o form an admissible 
basis of W. We can expand *$>w( x i z ) hi a formal series 

* w (x, z) = e xz ^ + E a *V) z ~* S j ■ ( L1 ) 

The wave function fy\v(x,z) can also be written in terms of the so-called wave 
operator Kyy- This is a pseudo-differential operator defined by 



K w {x,d x ) = l + Y j a k (x)d- k . (1.2) 



fe=i 

Then obviously 

#wr(a?,z) = ^(x,^)e xz . 
Introduce also the pseudo-differential operator 

P = K w od x oK-\ (1.3) 

For the treatment of the bispectral problem the following identity is crucial: 

P^w(x, z) = z^ w (x, z). (1.4) 

When it happens that some power of P, say P N , is a differential operator L we get 
that fyw( x i z ) is an eigenfunction of an ordinary differential operator: 

LV w (x, z) = z N ^ w (x, z). (1.5) 

It is easy to show that P N is differential iff the plane W is invariant under the 
multiplication by z N : 

z N W C W. (1.6) 

The submanifold of Gr consisting of planes satisfying ( |1.6| ) is denoted by Gr^ N \ 
A very important object connected to the plane W is the algebra A\y of poly- 
nomials f(z) that leave invariant: 

A w = {f(z)\f(z)W CW}. (1.7) 
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For each f(z) £ Aw one can show that there exists a unique differential operator 
Lf(x,d x ), the order of Lf being equal to the degree of /, such that 

L f (x, d x )^ w (x, z) = f(z)^ w (x, z). (1.8) 

Explicitly we have 

L f = K w of{d x )oK^\ (1.9) 
We denote the commutative algebra of these operators by Aw-, i-e. 

A w = {Lf\L f V w = f* w , f e A w }. (1.10) 

Obviously, Ay/ and Aw are isomorphic. We call Aw spectral algebra corresponding 



to the plane W . Following I. Krichever (see e.g. [KrN]), we introduce a rank of Aw 
to be the dimension of the space of eigenfunctions ^w\ this number is equal to the 
greatest common divisor of the orders of the operators Lf. At the end we define the 
spectral curve corresponding to the plane W to be Speo4jy (for definition see e.g. 
[ AMcDfl ). It is known that SpecA^ is an algebraic curve (see [|BC] , |, KrNfl ) . 



Remark 1.1 If ^>w{ x i z ) is we h defined for x = xq we set Vk = e x ° z z k and consider 
the subspace W x ° of V with an admissible basis = d kl $>w( x , z )\x=x - The wave 
functions of W x ° and W are connected by ^w x o(x, z) = e~ x ° z ^w(x + xq,z) and 
obviously 

K w *o (x, d x ) = K w {x + xq, d x ). 

These shifts are inessential for the bispectral problem and for our proofs. Through- 
out the paper we shall sometimes work with ^w x o calling it by abuse of notation a 
wave function of W and denoting it ^>w- 



Remark 1.2 In Sato's theory of KP hierarchy one usually considers the wave func- 
tion 



z~ k 



k=l 



depending on all times t\,t2, . ■ ■ (here t\ = x). Then the operators Kw(t,d x ) and 
P{t,d x ) are given again by formulae (|1.2| , |Q| ) with ak(x) replaced with afc(t) and P 
satisfies the following infinite system of non-linear differential equations 

where P k stands for the differential part of the fc-th power P k of the operator P, 
[. , .] is the standard commutator of pseudo-differential operators. The equations 
( |1.11| ) are called the KP hierarchy. If the plane W lies in Gr^ N ^ then ( |1.11| ) can be 
written in the form 

(L = P N ) and is called an N-th reduction of the KP-hierarchy or an N-th Gelfand- 



Dickey hierarchy 
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1.2 Darboux transformations 

In this subsection we recall the notion of Darboux transformations on objects con- 
nected to points of Sato's Grassmannian, introduced in our recent paper [ BHY3[ , 

Recall that a Darboux transformation [Da] of an ordinary differential operator L 
is obtained by presenting it as a product and exchanging the places of the factors: 

L = QP^L = PQ. 

A (monic) operator L is completely determined by its kernel: if {$0) • • • > l} is a 
basis of KerL then (see e.g. [|) 

_ Wr^p,...,^.!,*) 
L *" Wr^o,...,^) 

where Wr denotes the Wronski determinant. The next lemma answers the question 
when the factorization L = QP is possible (see e.g. ||). 

Lemma 1.3 (i) For a given basis {<&o, . . . , & n -i} of KerL set 

qx(x) = cUog$ , (1-14) 

Then the operator L can be factorized as follows 

L = (d x - q n )(d x - g n _i) ■ ■ ■ {d x - q x ). (1.16) 
(ii) L can be factorized as 

L = QP iff KerP C KerL. (1.17) 

In this case 

KerQ = P(KerL). (1.18) 

A slightly more general construction is the following one. For operators L and P 
such that the kernel of P is invariant under L, i.e. 

L(KerP) C KerP (1.19) 

we consider the transformation 

L^L = PLP-\ (1.20) 



The fact that L is a differential operator follows from Lemma |1.3| (ii). Indeed, 
L(KerP) C KerP is equivalent to KerP C Ker (PL). 

In JBHY3 ] we defined a version of Darboux transformation on points in Sato's 
Grassmannian and on related objects - wave functions, tau-functions and spectral 
algebras. 
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Definition 1.4 We say that a plane W (or the corresponding wave function *$>iy(x, z)) 
is a Darboux transformation of the plane V (respectively wave function ^fy(x, z)) iff 
there exist monic polynomials f(z), g(z) and differential operators P(x, d x ), Q(x, d x ) 
such that 

* w (x, z) = -^rP(x, d x )V v (x, z), (1.21) 
*v(x, z) = JJzjQ( x > d ^w(x, z). (1.22) 
An equivalent definition is that W is a Darboux transformation of V iff 

fV CW C-V (1.23) 

for some polynomials f(z), g{z). 

Simple consequences of Definition |1.4| are the identities 

PQ* w (x, z) = f(z)g(z)^ w (x, z), (1.24) 
QP^ v (x, z) = f(z)g(z)V v (x, z). (1.25) 

The operator L = PQ G Aw is a Darboux transformation of L = QP G Ay- 

Having in mind applications to the bispectral problem, the most important for 
our study is the case when 

A v = C[z N ], A V = C[L V ] (1.26) 

for some natural number N and a differential operator Ly of order N. (This is the 
simplest case of rank N spectral algebra with a spectral curve C.) Then due to 
(p.,25 ) we have 

f(z)g(z) = h(z N ), (1.27) 
QP = h{L v ) (1.28) 

for some polynomial h(z). In [|BHY3[| we connected the spectral algebra Ayy (re- 
spectively Aw) with Ay (respectively Ay). 

Proposition 1.5 (i) The Darboux transformations preserve the rank of the spectral 
algebras, i.e. if W is a Darboux transformation ofV then r&nkAw = rankAy. 
(ii) If Ay = C[L V ], ordLy = N then 

A w = {u G C[z N ] \ u{L v )KerP C KerP} (1.29) 

and 

A w = {Pu(L v )p- 1 | u G A w ) ■ (1.30) 
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1.3 Bessel operators, Bessel planes and related objects 

Now we define the planes of the Sato's Grassmannian on which we shall perform the 
Darboux transformations. For (3 G C N such that 

A n N(N - 1) , 

E^ = ~S — - ( L31 ) 

i=i 

we introduce the ordinary differential operator 

Pp{D z ) = (D z - fa)(D z - fa) ■ ■ ■ (D z - (3 N ), (1.32) 

where D z = zd z , and consider the differential equation 

Pp(D z )$p(z) = z N $>p{z). (1.33) 

For every sector S with a center at the irregular singular point z = oo and an angle 
less than 2tt the equation (1.33) has a solution <J?^ with an asymptotics 

oo 

Spiz) ~ Vp(z) =e z (l + Y, MP)z~ k ) (1-34) 

fc=l 

for \z\ — ► oo, z G S (see e.g. |Wa| ). Here a,k(/3) are symmetric polynomials in 
The function $3(2) can be taken to be (up to a rescaling) the Meijer's G-function 



*/»(*) = gS8((^) ^) (1.35) 



-2 



iV 



1 



TV 



- see |BE|, §5.3. 

The next definition is fundamental for the present paper. 

Definition 1.6 Bessel wave function is called the function ^fp(x,z) = ^fp(xz) (cf. 
HI |^]). The Bessel operator Lp is defined as 

L^x,^) = x~ N Pp{D x ). (1.36) 

A Bessel wave function ^> p defines a plane Vp G Gr (called Bessel plane) by the 
standard procedure: 

Vp = span{d^p(x,z)\ x=1 }. 

(In fact \l/y„(x, z) = e~ z ^ p(x+l, z), cf. Remark p~Tl]; we took xq = 1 because ^p(x, z) 
is singular at x = 0, arbitrary xo ^ will do.) 

Because ^(x,^) depends only on xz, i.e. 

=D^(x,if), (1.37) 

it gives the simplest solution to the bispectral problem: 

Lp(x, d x )Vp(x, z) = z N *p(x, z), (1.38) 
Lp(z, d z )Vp(x, z) = x N Vp(x, z). (1.39) 
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1.4 Involutions in Sato's Grassmannian 



In this subsection, following G. Wilson [Wl], we define several involutions on points 
of Sato's Grassmannian and on related objects - wave functions and wave operators. 



Besides the general properties of the involutions taken from [Wl], we specify their 
action on Bessel planes. 



Introduce after [DJKM] the non-degenerate form in V (realized as the space of 
formal Laurent series in z^ 1 ) 

B(f,g) = -Resoc f(z)g(-z)dz, f,g eY. 

If V 6 Gr is a plane, define aV £ Gr to be the plane orthogonal to V with respect 
to the form B( . , . ), to wit 



aV = {g(z) | B(f,g) = 0, for all f EV}. 



(1.40) 



Obviously, a(aV) = V, i.e. the map a is an involution. Following [Wl], we call it the 



adjoint involution. On the wave operator Kw fll.2| ) the involution a acts as [ Wl ] 



K aV = {K{ 



(1.41) 



where * is the formal conjugation on pseudo-differential operators, i.e. the antiau- 
tomorphism defined by d* = —d x , x* = x. For our purposes the most important 
property of the involution a is that it inverses inclusions, i.e.: 



if W C V, then aW D aV. 



(1.42) 



The following proposition will be used in the description of the action of a on Dar- 
boux transformations. Its simple proof is similar to that of Corollary 7.7 from [Wl]. 



Proposition 1.7 (i) If ^w(x, z) 



1 



■P(x,d x )^v(x, z), then 



^av(x,z) 



P*(x,d x )^ aW (x,z) 



(1.43) 



where g(z) = g(—z). 



(ii) flWlH Let Ay he the algebra of operators (|1.10|) . Then A a v consists of the 
conjugated operators of Ay . 

Proof. We have * w (x,z) = g- 1 (z)P^ v (x, z) = g' l {z)PK v e xz = PK v g~ 1 (d x )e xz 
which implies Ky/ = PKyg~ l (d x ). Applying the involution a we obtain 

K aW = {K* w y l = (P*)- 1 (K^)- l g(d x ), 



yielding P*K aW = K aV g(d x ) and hence (|L43|) . □ 

The sign involution s [ |W1[[ is defined on the wave functions by the property 



* s y(x,z) = V v (-x,-z). 



(1.44) 
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On the wave operators Ky(x,d x ) ( (1.44) ) translates into 

K sV (x,d x ) = K v (-x,-8 x ). (1.45) 

The plane sV is defined by 

sV = {f(-z) | f{z) G V}. (1.46) 

In the subsequent chapters we shall need the action of a and s on the Bessel 
planes. We describe them in the next proposition. 

Proposition 1.8 The involutions s and a act on Bessel planes Vp (/3 S C N ) as 
follows 

sVp = Vp, (1.47) 
oVp = V m , (1.48) 

where a{(5) = (N - 1)6 - 0, S = (1, 1, . . . , 1). 

Proof. We compute the action of the involutions on the corresponding wave func- 
tions. Obviously, 

^sV (x,z) = V p (-x,-z) = ®p(x,z), 

showing ( 1.47[) . 

Proposition L7, eq. ( 1.38 ) and the fact that = (—l) N L a ^ imply ^ a v^(x, z) = 
j{ z )^V a rm ( x i z ) f° r some formal power series 7(2;) in z. To show that 7(2) = 1 we 
notice that ^ a v (x,z) depends on xz. Indeed, ( 1.37j ) implies Ky^cx ^ c~ x d x ) = 
Kvp(x, d x ) for all c ^ and the same is true for Kav^ = (Ky^) -1 . □ 

We end this section by recalling the bispectral involution b which Wilson [Wl] 
introduced for the purpose of the bispectral problem. Contrary to a and s, the 
bispectral involution b is not defined on the entire Grassmannian. Whenever one 
can define b, put 

V bV {x,z) = Vv(z,x), (1.49) 

i.e. the involution b interchanges the places of the arguments x and z. 

A simple example of a point V € Gr where the involution b is well defined is the 
Bessel plane Vp. It immediately follows from the definition that 

<Sp(x, z) = ^p(z, x), i.e. bVp = Vp. 

In terms of the bispectral involution our approach to the bispectral problem can 
be formulated geometrically as follows: 

Find points V E Gr such that 

1) g{z)V C V for some nontrivial polynomial g(z); 

2) bV exists and f(z)bV C bV for some nontrivial polynomial f(z). 

A very important general property of 6, which we intend to use, is its connection 
to the other involutions [Wl]: 

ab = bas. (1.50) 
For completeness we also point out that the involution s commutes with a and b: 

as = sa, bs = sb. (1-51) 
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2 Polynomial Darboux transformations of 
Bessel wave functions 



The main purpose of this section is to introduce the submanifolds (denoted below by 
Gr^) on which, as we prove in the next section, the bispectral involution b is well- 
defined, and whose points correspond to bispectral operators. The points of Gr^ 
are obtained by a version of Backlund-Darboux transformation performed on Bessel 
wave functions (or equivalently, on polynomials h(Lp) of Bessel operators). Below 
we call these transformations polynomial Darboux transformations. Definition 



where this is done and the statement of Theorem 2.7, where we provide an equivalent 



definition, form the heart of the present section. Definition 2.5 has the advantage to 



be more natural and to supply an algorithmic procedure for constructing bispectral 



operators. The second definition (from Theorem 2/7) is more suitable for the proof 
of our bispectrality theorem in Sect. 3. 

The reader who wishes to see as soon as possible the main results of the paper 



can use the second definition, the one from the statement of Theorem 2.7, skipping 
its technical proof, which occupies half of the section. 

In the first half of this section we describe the kernel of the operator P from 
Definition 1.4 and (which is equivalent) - the conditions of the type as in [Wl], 



imposed on a Bessel plane, which define the corresponding Darboux transformation. 
To do so, we first need a description of the kernels of the operators h(Lp) which are 
polynomials of Lp. 

Fix (5 G satisfying ( |1.31 ) and let Vp be the corresponding Bessel plane (see 
Subsect. 1.3). Throughout this section W will be a Darboux transformation of Vp 



(we shall use the notation of Definition |1.4| with V = Vp). We shall need a lemma 
describing the kernel of the operator h(Lp) for an arbitrary polynomial h. 

Lemma 2.1 Let h{z) be a polynomial 

h(z) = z d ° (z - Xf) dl ■ ■ ■ (z - \?f r , Af ^ Af , A = 0, d, > 0. (2.1) 

Then we have 

(i) KeihiLfs) = 0[ =o Ker (Lp - \?) d \ 

(ii) (Lp) d = Lpd, where 

f3 d = ((3 1 ,f3 1 +N,...,(3 1 + (d-l)N,...,(3 N ,...,(3 N + (d-l)N). (2.2) 

(iii) If {Pi, . . . , /3n} = {ai, . . . , a%, . . . , a s , . . . , a s } with distinct oti,..., a s , then 



KevLp = span {^(ma;)*}^ 



l<i<s, 0<k<ki-l 

(iv) For A ^ 

Ker (Lp - X N f = span ^A.}^^ 
where e = e 2m l N is an N-th root of unity. 



14 



The proof being obvious is omitted (cf . Lemma |1.3| ) . □ 

Let us consider the simplest example of a Darboux transformation. Set 

h(z) = z d , g(z) = z n , f{z) = z dN ~ n (2.3) 

and 7 = (3 d , i.e. 

7(fc-l)d+j := Pk + U ~ l)N, 1 < k < N, 1 < j < d. (2.4) 

For an n-element subset I of {1, ... , dN} such that ji ^ 7-,- for i ^ j G I, we put 

KerP = span{x 7l } ig/ . (2.5) 

Such P corresponds to a Darboux transformation Vl/j(a;, z) of ^b{x, z). The following 
simple fact will be useful in the sequel. 

Lemma 2.2 ^j(x,z) is again a Bessel wave function: 

&l(x, z) = ^ 1+ dN8 I ^n8{x, Z). (2.6) 

Here and further we use the vectors 5j, 5 defined by 

(fr)i = (n V-Vr (2-7) 

] 0, if l $L I K ' 



and 



Proof. By definition 



1 for all i G {l,...,dN}. (2i 



¥j(a;, *) = z~ n L 7l * 7 (x, z), (2.9) 

where 7/ = {7i}j G /. Then ^j(x,z) is an eigenfunction of the differential operator 
L 7/ L 7 L~ , which is straightforwardly computed to be equal to L^^^-nS- 1=1 
Our next step is to study the spectral algebra Ap = Av of a Bessel plane Vp 
(Ol OCT)). 



sec 



Lemma 2.3 If L(x,d x )^p(x,z) = u(z)^ p(x, z) for some operator L € .4/3 and 
some polynomial u(z) € Ag ; i/ien L is a linear combination of Bessel operators L a , 
a eC k such that L a ^p(x,z) = z k ^lp(x,z). 

Proof. Let u(z) = Y^ u kZ k , Uk 7^ 0, < k < M for some M. Then for arbitrary 
c / we have 

u(cz)^ (x , z) = u(cz)^> (c _1 x , cz) = L(c^ 1 x,cd x ) y ip(x,z). 

This implies that u(cz) £ An and thus z k £ Ag. On the other hand D z Vp C Vp and 
the compatibility condition is of the form (D z — at\) ■ ■ ■ (D z — 0^)^(1, z) = z k ^> p(\, z) 
which implies L a ^ p(x, z) = z SS?g(x, z). □ 

Let us introduce the following terminology. We say that (3 € is generic if Vp 
is not a Darboux transformation of another Bessel plane Vw with (3' € C N , N' < N. 
The next proposition seems obvious but we do not know a simpler proof. 
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Proposition 2.4 For a generic (3 6 C we have 

Ap = C[z N ], Ap = C[Lp\. (2.10) 

Proof. We shall prove that if xankAp = r < N then there exists (3' € C (with 
Api = C[z r \) such that Vg is a Darboux transformation of Vw. The main idea is to 
apply to Vg Darboux transformations which lead again to Bessel planes but reduce 
the order of the operator Lp. Note that according to Proposition |l,5| they preserve 
the rank r of the spectral algebra. 

Split the set (3 into congruent modiVZ classes 




such that 

# - £ NZ for i + j and all s, t 

(N = N 1 + --- + N P ). 

By a Darboux transformation this (3 can be changed to 

iVi N p 

such that 

|Re(/?W - (3^)\ < N and ^ 0& for i ^ j (2.11) 

(see Lemma 

Suppose that ^ C[z N }. Then by Lemma [0] there exists a Bessel operator 
L a such that 

L a *fp(x,z) = z M y f3 ,(x,z), aeC M (2.12) 
and L a ^ L^, for any k. It is clear that 

L a L/3> = LprL a , (2-13) 

which is equivalent to 

{«! + N, a 2 + N, . . . , a M + N, , . . . , , . . . , . . . , } 

= {(3^ + M, . . . , (3^ + M, . . . , (3^ + M, . . . , [3^ + M,a u . . . ,a M } ■ 

Now if M > N this imply that f3' C a and therefore there exists a Bessel operator 
L a i such that 

L a = L a iLf)i and L a iLpi = LRiL a i. 

Repeating the same argument with a', we obtain that there exists L a satisfying 
(|2.13| ) with M < N. But then ( p.!3|) is equivalent to = V a . By Proposition 
|Lq r = rankAg = rankj4^/ = rankA Q divides M and N. If V a = Vp> = C[z r ] this 
finishes the proof. Otherwise we can repeat the above argument with V a instead of 
Vp. □ 
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Now we come to the main purpose of this section: the definition of manifolds of 
Darboux transformations, which will give solutions to the bispectral problem. To 
get some insight we shall consider, following Wilson [Wl], the geometrical meaning 



of Darboux transformations, provided by the so-called conditions C. 

Proposition [2.4| implies that for generic (3 eC N (|L27|, |L2|) hold with V = Vp 



and KerP is a subspace of Ker/i(L«). Each element / of KerP corresponds to a 
condition c (a linear functional on Vp), such that 

f(x) = (c,Vp(x,z)), (2.14) 

c acts on the variable z. These linear functionals form an n-dimensional linear 
space C {space of conditions) where n = ordP. In this terminology the definition of 
Darboux transformation can be reformulated as 

W = -^{v G V (c, v) = for all c G c\ 
g(z) I J 



(sec [Wl, BHY3| |). Following Wilson |W1[ |, we call the condition c supported at A 



iff it is of the form (cf. Lemma ^lj (iv)) 

%\z=X 



J2 akd 



'■' (2.15) 



(the sum is over k G Z>o and only a finite number of ^ 0). For Bessel wave 
functions this definition does not make sense when A = (since z) has a 

singularity at z = for N > 1). In this case we say that c is supported at z = iff 
it is of the form (cf. Lemma [2,1| (ii, iii)) 



The sums are over a G Ui=i{A + NZ>q} and < j < mult(a) — 1 where mult(a) is 
the multiplicity of a in the above union (only a finite number of b a j ^ 0). The space 
of conditions C is called homogeneous iff it has a basis of homogeneous conditions c 
(i.e. the support of c is a point). 

It is easy to see that if C is homogeneous then the spectral curve Spec^vi/ is 



rational and unicursal [Wl] (i.e. its singularities can be only cusps) - the condition 



c supported at A "makes" a cusp at A. For rank one algebras rationality and uni- 



cursality of SpecA^K are necessary and sufficient for bispectrality [Wl]. For rank 



N > 1 another necessary condition is that SpecAjy be T L^-invariant 1 i.e. 

A w C C[z N ]. (2.16) 

When W is a Darboux transformation of a Bessel plane Vg, with generic G C^, 
this condition is satisfied because of Propositions ^J|, |ll]. It is natural to demand 
that the space of conditions C (or equivalently KerP) also be Z^r-invariant. 
The Z 7v-invariance of KerP simply means that 

f{x) G KerP => f(ex) G KerP, e = e 2ni/N . (2.17) 
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It is easy to see that C is homogeneous and Z^v -invariant iff KerP has a basis which 
is a union of: 

(i) Several groups of elements supported at of the form: 

fc mult(ft+fc/V)-l 

9 1(J2 E b kj x^ +kN y j ) ] , 0</<io, (2.18) 



fc=o i=o 



where jo = max{j|frfcj ^ for some A;}; 

(ii) Several groups of elements supported at the points e l \ (0 < i < N— 1, A 7^ 0) 
of the form: 

fc 

J> fc ^^(^)L= £ *A, 0<i<iV-l. (2.19) 
fc=o 

Instead of ( |2.19| ) we can also take 

feo 

J2akD k z *p(x,z)\ z=E i X , 0<i<N-l, (2.20) 

k=0 

Denote by uq the number of conditions c supported at (i.e. the number of elements 
of the form (2.18) in the above basis of KerP). For 1 < j < r denote by nj the 
number of conditions c supported at each of the points £ l Xj, < i < N — 1 (i.e. the 
number of groups of elements of the form (2.19) with A = Xj). 
We have at last arrived at our fundamental definition. 

Definition 2.5 We say that the wave function $?w(x,z) is a polynomial Darboux 
transformation of the Bessel wave function ^^(x,z), (3 6 C , iff (|1 .21 ) holds (for 
V = Vp) with P(x, d x ) and g{z) satisfying: 

(i) The corresponding space of conditions C is homogeneous and Z7v-invariant, 
or equivalently KerP has a basis of the form (2.18, 2.19j ). 

(ii) The polynomial g{z) is given by 

g{z) = z no {z N - Af ) ni ■■■(z N - \?) nr (2.21) 



where nj are the numbers defined above. 



We denote the set of planes W satisfying (i), (ii) by Grs(P) and put Gr B 
\J p Gr B (0), /? G C^-generic. 

We point out that the form ( |2.2l| ) of g(z) was introduced for N = 1 by Wilson 
[|Wll1 . (Note that g(z) = z n ° n^illilo 1 ^ ~ £%■)"'•) We make this normalization 
in order that ^bw( x , z ) = ^wi z -,x) be a wave function; for the bispectral problem 
it is inessential. 

Definition 2.6 We say that the polynomial Darboux transformation *&yv(x, z) of 
^>p{x,z) is monomial iff 

g(z) = z n ° 

(i.e. iff all conditions c are supported at 0). Denote the set of the corresponding 
planes W by GrMB(P) and put Gr^ B = U/3 Gr (/3) , £ C^-generic. 
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The next theorem provides another equivalent definition of Grs((3) and is used 
essentially in the proof of the bispectrality in the next section. 

Theorem 2.7 The wave function ^w(x, z) is a polynomial Darboux transformation 
of the Bessel wave function ^fp(x,z), (3 G C^, iff ( |L2g , [L22| \L2l\ , [O^) hold (for 
V = Vp) and 

(i) The operator P has the form 



p(x,d x ) = x~ n Y / Pk(x N )(xd x ) k , 



(2.22) 



k=0 



where pk are rational functions, p n = 1 . 
(ii) There exists the formal limit 



lim e xz ^ w (x,z) = 1. 



(2.23) 



The proof will be split into three lemmas. Before giving it we shall make a few 
comments. 

The rationality of P is always necessary for bispectrality [ DG , Wl |, Q2.22D 
also imposes the ZAr-invariance. The condition ( 2.23| ) is necessary in order that 
^bw( x > z ) = ^w( z , x ) be a wave function. The limit in (|2.23 ) is formal in the sense 
that it is taken in the coefficient at any power of z in the formal expansion ( |1 . 1| ) 
separately, i.e. 

lim aAx) = for all j > 1. (2.24) 

x— >oo 



Our first lemma is similar to Proposition 5.1 ((i) =4* (ii)) from | W1 |. 

Lemma 2.8 If P has rational coefficients and is 7L^ -invariant (see ( |2.22| )) then the 
conditions C are homogeneous and Z^r -invariant (see ( 2.18 , 2.19 )) 

Proof. If KerP = span{/o, . . . , f n -i}, the second coefficient of P is 

-a 2 logWr(/ ,...,/ n _i) 

and is rational. Lemma implies that Wr(/o 5 • • • ; fn-i) is of the form 

x a e Xx x (Laurent series in x _1 ). 

In particular each element of KerP is a sum of terms of the form 

e Xx x (Laurent series in x -1 ) or x a (lnx) k . 

We order the (finite) set of all such e Xx and x a (lnx) k occuring in KerP. The highest 
term in Wr(/j) is just the Wronskian of the highest terms of the /j. If it vanishes 
then the highest terms of the fi are linearly dependent, so by a linear combination 
we can obtain a new basis with lower highest terms. So we can suppose that the 
highest term of Wr(/j) is non-zero. Repeating the same argument with the lowest 
term, we shall finally obtain a basis whose elements consist of only one term, i.e. are 



homogeneous (cf. |Wl| ). 
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Because the coefficients of P are rational, implies that it does not matter 

which branch of the functions x a (lnx) k in KerP we take for x G C. Let 

30 
3=0 

with fj(x) = b a jX a . Then ^ fj{x)(\nx + 21t\i) 3 G KerP for arbitrary I £ Z and 
also for / G C since it is polynomial in /. Taking the derivative with respect to I we 
obtain that 

jo 

^2f 3 (x)j(lnxy- 1 

3=0 

also belongs to KerP. 

On the other hand the ZAr-invariance of P (see (|2.17 )) implies Yl fj( £x )(^ nx + 
2m/ N) 3 ' G KerP and 

jo 

^2 fj (ex) (In x) j G KerP 

3=0 

for e = e 27ri / N . 

Now it is obvious that KerP has a basis of the form ( [2.18| , 2.19| ). □ 

Lemma 2.9 If KerP has a basis of the form ( 2.1 8| , 2.19| ) then P has rational coef- 
ficients and is Z^v -invariant (see ( 2.22; )). 

Proof. Consider first the case when the basis of KerP is 

fi(x)=Y,°k%*f}(jx,z)\ t=x , 0<i<N-l, A^O. 

k 

We shall show that det (d™ 3 fi(x)) 0<i j <N _ 1 is a rational function of x for arbitrary 
nj G Z> . Using (|L37j |T3|) we can express all derivatives of ^(a;,^) (both with 
respect to z and x) only by ^^(a;, z), < A; < N — 1, to obtain 

iV-l 

'./)!■<•) = ^ « fci (x, AR fc f ^x, A) (2.25) 

fe=0 

with rational coefficients ay. Therefore 

det $?fi(x)) = det(a kj (x, A)) det^^x, A)). 

But det (5^ ^/3(e l x, A)) = const because the second coefficient of Lp — A is 0. If the 
basis /o, . . . , fmN-i of KerP contains m groups of the type considered above (i.e. 
(|2.19D ) we can represent the matrix 



(^7i(aO) <i,,-<miV-l> %GZ>o, 
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in the block-diagonal form 





V 





w 2 



W m J 



where each block W s has the form already considered above. This can be achieved 
by columns and rows operations, using the representation ( |2.25| ). 

If in addition there are some groups of elements of the form (|2.18j ), we kill the 
logarithms by columns operations and then cancel the powers ar* from the numerator 
and the denominator of ( |1.13| ). □ 

Lemma 2.10 If C is homogeneous and h, g are as in (|2.1|), ( |2.21| ) ; then ( |2.23| ) is 
satisfied. Conversely, ( 2.23|) implies ( 2.21|) . 

Proof. The second part of the lemma is an obvious consequence of the first one. 

For a basis {<&i(x)}o<i<dN-i of Ker/i(L^) (d = degh) we consider the basis of 
KerP 



dN-l 



fk(x) -- 

Formulae ( [1.21 , 1.13 ) imply 



^2 a hi <&i(x), < k < n-1. 



i=0 



Wr(/o(x), ■ ■ ■ , / w _i(g), *f}(x, z)) 

0(z)TO(/ o (x),...,/n-l(aO) 
^det^ 7 Wr(^(x))^/(j;,z) 



(2.26) 

(2.27) 
(2.28) 



£det^Wr($/(x)) 

The sum is taken over all n-element subsets 

I = {i <h < ... < i n _i} C {0, 1, . . . , dN - 1} 

and here and further we use the following notation: A is the matrix from ( 2.26| ) and 
A 1 = (a kt i t ) <k,i<n-i is the corresponding minor of A, $/(x) = {® io (x), . . . , ^^(x)} 
is the corresponding subset of the basis {&i(x)} of Ker/i(L / g) and 



Wr($j(g),^(a;,z)) 
g(z)Wv{^ I (x)) 



(2.29) 



is a Darboux transformation of \P p{x, z) with a basis of KerP fk = 

Using ( 2,28j ) it is sufficient to prove (|2.23 ) for ^i(x, z), hence we can take KerP 
consisting of functions 



f i (x)=d k z ^ f3 (x,z)\ z=Xi , 0<i<p 
fi(x) = x ai (\nx) li , p < i < n — 1. 



(2.30) 
(2.31) 



We shall consider the case when Aj ^ \j for i ^ j. The general case can be reduced 
to this by taking a limit. In the formula ( |2.27 ) we expand the determinants in the 
last n — p columns (using the Laplace rule) : 



Wi(f,*p) = Y ( ±deb{di>f i (x),dk9f i (x,z)) o< s < P -det (dt hix)) P+1 < s < n ; (2.32) 

*■ ' 0<i<p-l p<i<n~l 
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Wr(/) = £±det (^/i(*))o<.,i<p-i-dflt (^/i(^)) p < s , i < n _ 1 , (2-33) 

where the sums are over the permutations (jo, ... , j n ) (resp. (jo, . . . , jn-i)) of (0, ... , n) 
(resp. (0, . . . , n - 1)) such that j < . . . < j p and j p+ \ < . . . < j n (resp. jo < ■ ■ ■ < 
jp-i and j p < . . . < j n -i). We extract the terms with the highest power of x in the 
numerator and in the denominator of ( |2.27 ). Obviously, 



det {di a fi{x)) p<i<n l = const • x^=p ai ^ sja Rj(lnx) (2.34) 

for some polynomials Rj(\nx) (J is the permutation (j s )). On the other hand for 
< i < p- 1 

dtfi(x) = dl°d k z ^p(x,z)\ z=Xi = x k >e x > x (Af + 0{x- v )) (2.35) 

and 

dfrVpix, z) = e xz (z js + Oix- 1 )) . (2.36) 
Now it is easy to see that the leading terms are obtained for the permutations 

(n — p, n — p + 1, . . . , n, 0, 1, . . . ,n — p — 1), 

respectively 

(n — p, n — p + 1, . . . , n — 1, 0, 1, . . . , n — p — 1). 



Substituting (fO|, g]3§ fO§) in p^2] g3g ) and canceling the determinant ( gjgg ) 



for J = (0, 1, . . . , n — p — 1), we derive that 

lim e~ xz P(x,d x )ty l3 (x,z) 

x^oo 

is a fraction of two van der Monde determinants and therefore is equal to g(z). □ 

3 Bispectrality of polynomial Darboux transformations 



In this section we prove the main result of the paper, Theorem |3.3| , claiming that 
polynomial Darboux transformations (see Definition [2.5| ), performed on Bessel oper- 
ators, produce bispectral operators. On its hand Theorem fO is an almost obvious 



consequence of Theorem 3.2 in which we prove that the bispectral involution is 



well-defined on the submanifolds Gr^(/3) and maps them into themselves. The im- 



portance of Theorem 3.2 is not only to provide a proof of our main result (Theorem 
but also to enlighten the bispectral involution. Its proof uses only the definition 
of polynomial Darboux transformation from Theorem [O] (i.e. it does not use Def- 
inition ^5| ). On the other hand, the proof is completely constructive and together 
with Definition |2.5| it provides an algorithmic procedure to compute bispectral wave 
functions and the corresponding bispectral operators. This procedure is described at 
the end of the section. Many examples computed by making use of it are presented 
in Sect. 5. 
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Let Vg be a Bessel plane for a generic P G (i.e. Vp is not a Darboux trans- 
formation of Vpi with P' G , N' < N). In this section W will be a polynomial 
Darboux transformation of Vp, i.e. 

W e Gr B ((3). 

We use the notation from ( 1.21| , 1 .22) ) with V = Vp. In the next proposition we 
show that the manifold of polynomial Darboux transformations is preserved by the 
involutions a and s (introduced in Subsect. 1.4). 

Proposition 3.1 IfWe Gr B (P), then 

(i) sW e Gr B (p); 

(ii) aW G Gr B {a{p)), where a((3) = (N - 1)6 - P, 6 = (1, 1, . . . , 1). 



Proof. First recall that (Proposition |1.8[ ) sVp = Vp and aVp = V a m\. We shall study 
the action of the involutions on ^f\y(x, z) and check that the conditions of Theorem 
^7?] are satisfied. 

(i) is trivial because 

^sw(x, z) = V w (-x, -z) = l _ , P(-x, -d x )^p(x, z). 

9\ z ) 



To prove (ii) we note that the Z^y-homogeneity of P (see (2.22)) is equivalent to 

P(ex,e- 1 d x )=e- n P(x,d x ), (3.1) 

for n = ordP, s = e 2m / N . It follows from ( |1.28| ) that the operator Q (from ( 1.22[ )) has 
the same property and also that Q = h(Lp)P~ 1 has rational coefficients. Proposition 
|0] implies that ^* a w is a Darboux transformation of ^ a (/3) with 



1 

9~z) 



^ aW (x,z) = Tr - x Q*(x 1 d x )^ a ( K p ) {x,z). 



Obviously, Q* also satisfies Q3.1| ), To check ( |2.23j ), we set 

oo 

K w = 1 +^2,aj{x)d~ 3 

3=1 

(see Q ED)- Recalling that 

oo 

^W = l + £(-^)~ i «i(a?) 

3=1 

and 

oo 

K aW = 1 + £ V = W)" 1 , 

i=i 

we compute the coefficients bj(x) inductively and find that all of them are poly- 
nomials in aj(x) and their derivatives. But by Theorem 2/7 all a,j(x) are rational 



23 



functions of x and lim. E _ >00 a,j{x) = 0, which leads to hin^oo bAx) = for all j > 1. 
This proves ( gjg) for aiy (cf. □ 

Proposition 3.1 shows that the involutions a and s preserve Gr^\ The central 
result of the present paper is that the bispectral involution b has the same property. 
It immediately implies that wave functions *$>w with W G Gr% give solutions to 
the bispectral problem. Our next theorem addresses this issue. 

Theorem 3.2 IfW G Gr B ((3) then bW exists and bW G Gr B (f3). 

Proof. Before proving the existence of bW, we shall find an analog of ( |1.21| ) for 
^bw(x, z) = ^>w{z, x), i.e. we shall show the existence of an operator Pb{x, d x ) and 
a polynomial g\> (z) such that 



9b(z) 



P b {x,d x )^p{x,z). 



From ( p.22 ) it follows that the operator P can be written as 

1 n 

P(x '^ ) = ^7^T^^ (x7V)( ^ )fe ' 

where now p k (x N ) are polynomials. Use (L37-|L3^) to obtain 

1 



(3.2) 



(3.3) 



^ w (x,z) 



x n p n (x N )g(z) 
1 

x n p n {x N )g{z) 



Y,Pk{x N ){xd x ) k *p{x,z) 
J2(zd z ) k p k (L (z,d z ))^p( 



x,z) 



This implies ( |3.2[) with 



1 n 

Pb(x,d x ) = — - y2(xd x ) k p k {Lp(x,d x )), 



k=0 



(3.4) 
(3.5) 



g h (z) = z n p n (z N ). 

Now we can prove the existence of bW, i.e. that ^f b \y(x, z) is a wave function 



(see (|1.1|) ). Indeed, using ( |3.2j ) we can differentiate the formal expansion (1.34) of 
^p(x,z) = ^p(xz); expanding g^ 1 (z) at z = oo we obtain 

■$ bW (x,z) = e xz h{x)z- k 

k>ko 

for some finite k$. Note that the coefficients b k {x) are rational. On the other hand 



xzx aAz)x~ 3 



^ bW {x, z) = ^w(z, x) = e xz ^ WJ 

i>o 



with rational aj(z) such that (see ( gjg )) 

lim Oj(z) = 0, j > 1; oq(z) = 1. 



(3.6) 
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These two (formal) expansions of ^hw( x , z) are connected by 

aj(z) = ^2 b kj z~ k , 
k>ko 

^2b k jX~ j , b kj = for j < 0. 

3 



b k (x) 



where 

Now fl3.6[ ) implies 
This shows that 



is a wave function. It is clear that it satisfies (|2.23| ) as well. 
To show an analog of ( |1.22| ), i.e. that 

1 



b k j = for k < j > 1. 
^ bW {x,z) =e xz (l + ^b k {x)< 



k>l 



fb(z) 



Qb{x,d x )^ bW (x,z) 



(3.7) 



with an operator Q\> and a polynomial /b, we shall use the above proven identity 
;ead of W . I 

^ asW (x,z) 



•2J) with asW instead of W . It follows from Proposition |L7j that 

1 



Q*{-x,-d x )^ a{p) {x,z). 



(3.8) 



Proposition 3.1 and Theorem |2.7| (i) allow us to present Q*(—x, —d x ) in the form 

(3.9) 



Q*(-x, -d x ) = -^^n\ J2 Qs(x N )(xdx 
with polynomials q s (x N ). Then 



s=0 



^basw(x,z) 



f(x)z™q m (z N ) 



J2(xd x ) s q s (L am (x,d x ))V 



a(/3) 



,X,Z) 



(3.10) 



The identity ab = bas [ |W1[| and Proposition L7 now lead to ( |3.7D with 

Q h (x,d x ) = 

= f2Qs((-l) N Mx,d x ))(-xd x -lYj^ 



jj-j ^2(xd x ) s q s (L am (x, d x )) 



s=0 



s=0 



and 



(3.11) 

J {X) 

f h {z) = {-z) m q m ((-z) N ) . (3.12) 

From (2.21) and ( [3.4| ) it is obvious that is Z^y-homogeneous. This completes the 
proof of Theorem [T^. □ 

An immediate corollary is the following result, which we state as a theorem 
because of its fundamental character. 
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Theorem 3.3 If W € Gr^ V) then the wave function ^fw(x,z) solves the bispectral 
problem, i.e. there exist operators L(x,d x ) and A(z,d z ) such that 



L(x,d x )V w (x,z) = h(z N )^ w (x,z) 
A(z,d z )* w (x,z) = @(x N )V w (x,z) 



Moreover, 



rankAw = rankA^iv = N. 

Proof. (|3~13[ , follow from flPIj HHI U, U) if we set 

L(x, 9a;) = P(x, d x )Q(x, d x ), h{z N ) = f(z)g(z); 
A(z, d z ) = P h (z, d z )Q h (z, d z ), Q(x N ) = f h (x)g h (x) 

The eq. ( |3. 15| ) follows from Propositions |1.5| (i) and [2.4[ □ 



(3.13) 
(3.14) 

(3.15) 



(3.16) 
(3.17) 



Example 3.4 All bispectral algebras of rank 1 are polynomial Darboux transfor- 
mations of the plane H + 
Bessel with 



{z k }k>o ( see |W1| |). This corresponds to the N = 1 



/?=(0), L (0) =d x , V {0) = H + = {z k } k > , ^ (0) {x,z) = e xz . 
Every linear functional on H + is a linear combination of 

e(k,X)=d k z \ z=x 

and h (L(p)) = h(d x ) is an operator with constant coefficients. The "adelic Grass- 

mannian" Gr ad , introduced by Wilson |W1], coincides with GV#((0)) (= Gr^). In 
our terminology the result of | W1 | can be reformulated as follows. 

All bispectral operators belonging to rank one bispectral algebras are polynomial 
Darboux transformations of operators with constant coefficients. □ 



Remark 3.5 The eigenfunction tyw(x,z) from eq. ( [1.21 ) is a formal series. Let 
$p(x, z) = $>p(xz), where &b(z) is the Meijer's G-function ( 1.35| ) (or any convergent 
solution of ( |l .33|) in arbitrary domain) and set 



Then 



®w(x,z) 



1 



g{z) 
i 



9b {x 



P(x,d x )$i3(x,z). 

-P h (z,d z )^f3(x,z) 



(3.18) 



(3.19) 



because of ( |1.33| ) and xd x <S}p(x, z) = zd z $>p(x, z). The equations QP = h(Lp) and 
QbPb = @{L /3 ) imply 



$/3(x,z) 
®p(x,z) 



1 



1 



fb(x 



■Q(x,d x )$p(x,z), 
■Qb(z,d z )$ w (x,z) 



(3.20) 
(3.21) 
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So, we proved that $wi x > z ) is a convergent bispectral eigenfunction of the same 
operators L(x,d x ) and A(z, d z ) as ^w{x,z). The involutions a, s and b can be 
defined on the manifold of "convergent" polynomial Darboux transformations ( |3.18| ) 
by the equations (|1.43 , 1.44 , 1.4S| ) in which ^ is replaced by $ and they preserve it 
(Proposition |1.7| (i) now becomes a definition). The validity of the equation ab = bas 
in the "convergent" case is a consequence of that in the "formal" one (see the proof 
of Theorem |3.2[) . The rationality of the coefficients of the operator P{x,d x ) implies 
that its kernel has one and the same form (see eqs. ( p. 18 , 2.19| )) in ^- and in $-bases. 
□ 



It is not difficult to provide an explicit algorithm for producing bispectral pairs 
L(x,d x ), A(z,d z ). Although obvious we have collected the steps of this algorithm 
as they are scattered in the present and the previous sections. 

Step 1. Choose an arbitrary set of conditions based in some points Ao = 0, Ai, . . . , A r 
of the form (|2T8| , i.e. a basis of KerP. The proof of Lemma |2.9| provides an 

explicit computation of the coefficients of P in terms of the coefficients a&, bkj in 
KerP. The polynomial g(z) is given by Definition |2.5| (ii). 

d. 



Step 2. Take h(z 



y d N 



n 



i=i 



A 



N 



with high enough powers do,. . . ,d r 



such that KerP C Ker/i(L^) (cf. Lemma 2.1). The minimal such d/s can be com- 
puted as follows. 

(i) For a condition, supported at 0, of the form ( 2.18Q set j(k) = max{j\bkj ^ 0}, 
< k < k . Let Pi + kN = (3 is + p s N for < s < mult (A + kN) - 1 with 
<Po < ■■■ < Pmuit(ft+fcAf)-i and i s / i t for s^t. Then set 

d = 1 + maxp j(lt ), 



the maximum is over all k and all conditions of the form (2.18). 

(ii) For a condition, supported at Xj ^ 0, of the form ( 2.19 ) let ko = max{A;|afc ^ 
0}. Then set 

dj = 1 + max ko , 



the maximum is over all conditions of the form ( 2,19| ) supported at Xj. 
Then put f(z) = h{z N ) / g(z). 

Step 3. Find the coefficients of the operator Q(x, d x ) recursively out of the equation 
Q(x, d x )P(x, d x ) = h{Lp{x,d x )). Then L(x,d x ) = P(x,d x )Q(x,d x ). A lower order 
operator L can be constructed using Proposition |l.5| , i.e. find u(Lg) such that KerP 
is invariant under u{Lp) and then L out of the equation LP = Pu(Lp). 

Step 4. Compute by Q P h (x,d x ) and by Q g h (z). Also fml) and (|3T2| ) give 
Qb{x,d x ) and fh(z). All expressions are explicit in terms of the coefficients of the 
operators P and Q. Then A(z,d z ) = Pb(z, d z )Q^(z, d z ) and ©(x) = fb(x)g\ ) (x). 



4 Polynomial Darboux transformations of Airy planes 



This section contains analogs of the results from Sections 2 and 3 but here the 
building blocks are (generalized) Airy operators (see [KS, pij| ) instead of Bessel 
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ones. There is a minor difference in the organization of the present section compared 
to that of Sections 2 and 3. Here we give the definition of polynomial Darboux 
transformations on Airy wave functions (see Definitions 4.2, 4.3) in the spirit of the 
one provided by Theorem 2.7 . Then we prove our main result Theorem [4.5| (which 
is an analog of Theorem 3.2). As in Sect. 2, it automatically implies bispectrality of 
the polynomial Darboux transformations. At the end, in Proposition |4.9| we show 
that Definition 



is equivalent to a second one (analog of Definition 2.5) in terms 
of conditions on Airy planes. This is again important for algorithmic computations, 
some of which are presented in the next section. 

First we recall the definition of (generalized higher) Airy functions. 
For a = (qq, oii, 03, • • • > QW-i) € C^ -1 , «o 7^ 0, consider the Airy operator 



N-l 



L a (x, d x ) =d£ - a x + ^ aid£ 1 = P a >{d x 



OtQX 



i=2 



where a' = (02,03, ■ ■ ■ , ajv-i)- The Airy equation is 

L a (x,d x )$(x) = 0, i.e. P a ,(d x )$(x) = a x<S>(x). 



(4.1) 



(4.2) 



Example 4.1 When ckq = 1, a' = eq. ( |4.2[ ) becomes the classical higher Airy 
equation (cf. [KS]) 

d£ $(x) = x$(x). (4.3) 

In every sector S with a center at x = 00 and an angle less than Nir/(N + 1), it has 
a solution with an asymptotics of the form (see e.g. [|Wa[| ) 

iV+l 00 

$(x) ~ x~ i ^re^i x ^ T ~ (l + ^ajX _i//Ar ^, \x\ -> 00, x G 5. (4.4) 

i=l 

□ 



Similarly, in each sector S 1 as in Example [O] eq. ([4^) has a solution with an asymp- 
totics of the form 



$(x) ~ * a (x) := x d / N e Q{xl/N) (l + Yl a i x ' i/N ) > M ^ oc, x e S (4.5) 

i=l 



for some d £ C and a polynomial Q(x) of degree N + 1 with leading coefficient 
Hojj^iX N+1 , where ceo = Ho ■ The solution $ is by no means unique, but d, Q and 
all cii are uniquely determined and do not depend on S. In the sequel we shall deal 
only with ty a , which is a formal solution of eq. (|4,2j). 



Definition 4.2 For each a € 1 we call an Airy wave function the following 
function 

V^(x, z) := 4z~ d e- Q ^ lz ^ a (x, z), (4.6) 

where 

# a (x, z) := ^ a {a.Q 1 z N + x). 
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It is easy to see that ip a is indeed a wave function if we expand ^> a (a z + x) at 
x = 0: 

(a^z N + x)-^ = J] (~ l/ k N ) {^ l zr~ kN x k (4.7) 
fc>0 ^ ' 

(we shall always use \xq as an iV-th root of ao). 

The plane in Sato's Grassmannian corresponding to ifj a (x,z) will be called an 
Airy plane and will be denoted by V a . 

Obviously, *ff a (x, z) solves the bispectral problem 

L a {x, d x )V a (x, z) = z N V a (x, z) (4.8) 

L a (aQ 1 z N , d a -i zN )^ a (x, z) = a x^ a (x, z) (4.9) 

because 

d x V a (x, z) = d a -i zN V a (x, z). (4.10) 

It is clear that ifj a satisfies Q4.8|) and analogs of ( |4.9| , |4.10| ) obtained by conjugating 
by z~ d e~Q^o z ). (Up to this conjugation ^ a and i/) a give one and the same solution 
to the bispectral problem.) 

We shall define polynomial Darboux transformations of Airy planes as in the 
Bessel case (see Definition |2.5| and Theorem 2.7). Before that we shall define a bis- 
pectral involution b\ on them. Note that the involution b from |W1| (see Subsect. 1.4) 
is not well defined on V a (i.e. ip a (z,x) is not a wave function). The properties of b 
we would like b\ to have, are: 

1) it has to interchange the roles of x and z; 

2) it has to preserve Airy planes. 
Therefore we define 

biV a (x, z) := V a (x, z) = * a (a 1 z n , fi x 1/N ), (4.11) 

or equivalently, 

hijj a (x, z) := ijj a {x, z) = 4x d/N z~ d e Q ^ xl/N) - Q ^ lz ^ a {a a l z N ,^x 1/N ). (4.12) 

For a Darboux transformation W of V a we define ipbiW an d ^biW i n a similar way. 
(We still do not know whether b\W € Gr, the notation ip^w is still formal.) 

Definition 4.3 A Darboux transformation W of an Airy plane V a is called polyno- 
mial iff (in the notation of Definition |L4|) 

(i) the operator P has rational coefficients; 

(ii) g(z) =g 1 (z N ),g l eC[z\; 

(iii) lim e~ xz ijj blW (x, z) = 1. 

(The limit is formal and has the same meaning as in ( 2.23[) .) 

Denote the set of all such W £ Gr by Gr^(a) and put Gr^ = UaeC^- 1 ^rA(a). 
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Remark 4.4 The parts (i) and (ii) of the above definition remain the same if we 
substitute -0a and Vw by ^f a and ^w, where 

i(j w {x, z) := i4z~ A e- Q ^ lz) *w(.x, z). (4.13) 

The main result of this section is that Gr^(a) is preserved by the involution b\. 

Theorem 4.5 (i) IfW € Grj\.(at) } thenipb t w{x,z) is a wave function corresponding 
to a plane b\W £ Gr^(a). 

(ii) For a € C^ -1 the spectral algebra Ay a is C[L a ]. 

An immediate corollary is that the planes W € GrA{ot) give solutions to the bis- 
pectral problem of rankiV: 

rankAvK = vankA^w = N. 



The proof of Theorem is completely parallel to that of Theorem |3.2| , We shall 
be very brief, indicating only the major differences and the most important steps. 
We start with a lemma illuminating the purpose of the constraints (i) and (ii) in 
Definition (cf. Q). 



Lemma 4.6 IfW G Gr^(a), then 

V bl w(x, z) = —j-^Pbix, d x )V a (x, z), (4.14) 

gb { z) 

Pb is with rational coefficients and gb is polynomial in z N . 
Proof. We compute 

P(a^z N , d a ^ zN )^ a (a^z N , w 1/N ) ! 
= —, Tim = ^rPb(x,d x )^ a (x,z), 

where if 

n 

P{x,d x ) = — _^p fc ( x )^, g{z)= gi {z N ) (4.15) 

with polynomials pk and g\, then (using (4.8, |4.10D ) 

1 ™ 

P h (x, d x ) = — V d k xPk (a^L a (x, d x )), (4.16) 

giM ^ 

g h (z)=p n {a^z N ). (4.17) 

□ 

The proof that ipbiwi x > z) is a wave function is the same as in the Bessel case, 
using the above lemma and the condition (iii) of Definition 4.3. Now the identity 
ab = bas [Wl] is modified in the following way. 
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Introduce the maps p and p 1 as follows 

% w (x,z) := -^w^x 1 * ,hqZ 1/n ), 
■$ p -i w (x,z) := ^ w (p x 1/N ,ao 1 z N ). 

The notation pW, p~ 1 W is formal - these are not planes in Gr. But 

V blW (x,z) = ^ bpW (x,z) = V p -i bW (x,z) 

corresponds to the wave function ipfawfei z ) an d to b\W G Gr. 
Multiplying the identity ab = has on the right by p, we obtain 

ob\ = b\a\, where a\=p~ l asp. (4-18) 

Note that for W G Gta{oi) clW, b\W and hence a\W are planes in Gr. 

The next lemma gives the action of the involutions on the Airy planes (the proof 
is the same as that of Proposition |1.8|) . 



Lemma 4.7 (i) sV a = V a ( a \, where s(a) = ((— l) N+1 ao, ct2,— 03, . . . , (— 1)^ l a.N-i); 
(ii) aV a = aiV a = V a ^, where a(a) = ((-l) N a ,a 2 , -a 3 , . . . , (-l)^ -1 ^-!)- 

We also need an analog of Proposition jjP| . 

Lemma 4.8 IfW€ GrA{a), then aW and a\W belong to Gr^(a(a)). 



For the proof we need an analog of Proposition 1.7 for a\. A simple computation 
shows that if 

^v(x,z) = —-Q{x,d x )^ w {x,z) 

for V, W G Gr A {a) and 

Q{x,d x )=Y 4 1kix)^ 

then ^ 

^ aiW (x,z) = j^—Q* 1 (x,d x )^ aiV (x,z) 

with 

q* 1 = + ^s- 1 ) Vi) ( "- 1)fc ?*((-i)"*). 

The rest of the proof is left to the reader. □ 



The proof of part (i) of Theorem 4J3 is completed exactly as in the Bessel case. 
For the part (ii), we note that while the Bessel wave functions are 11 multiplication 
invariant" , the Airy ones are "translation invariant" . More precisely, for arbitrary 

c G C 

^ a (x + c, (z N -a c) 1/N ) = ^ a {a^ 1 {z N -a c) + x + c) = ^ a (a Q l z N + x) = $> a (x,z) 

(expand (z N -a c) l ' N = £ fc > (*f )z 1 - kN (-a c) k ). Let u(z) G A a , L(x, d x ) G A a 
and 

L(x, d x )^ a (x, z) = u(z)^ a (x, z) 
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(this is equivalent to Lip a (x,z) = uip a (x, z)). Then 

L(x + c, d x )V a (x, z) = u((z N - aoc) 1 ^)^, z) 

and L(x + c, d x ) G A a , u((z N — aoc) l / N ) G A a . But A a C C[z], therefore «((; 
aoc) 1 ^) G C[z] for all c and u(z) G C[z N ]. 



N 



This completes the proof of Theorem 4.5. 

At the end of this section we note that an equivalent definition of Gr^(a) can 
be given in terms of conditions C (cf. Sect. 2). 

Using the translation invariance of ^> a we can suppose that none of the conditions 
C is supported at 0. Then we have an analog of Theorem 2.7. 

Proposition 4.9 The Darboux transformation WofV a is polynomial iff 

(i) The space of conditions C is homogeneous and Z^v -invariant. Equivalently, 
KerP has a basis of the form 



fij( x ) = y^ j akid^Tp a (x,e j z)\ z= 
k 

= J2a ki e- jk d^ a (x,z) 

k 

< j < N — 1, 1 < i < r (for some r), Aj ^ 0. 

(ii) The polynomial g(z) has the form ( 2.2l| ), i.e. 



(4.19) 



g(z) = (z N -Af) ni ■■■(*" -A?) 

where n« is the number of conditions C supported at each of the points e-'Aj, < 
j<N-l. 

The proof of the "if" part is the same as in the Bessel case and will be omitted. (In 
fact, most of the proofs in Sect. 2 are valid in a more general situation.) 

The "only if" part is also similar to the corresponding result in the Bessel case 
but some more explanation is needed. For fixed A ^ we shall use representations 
of the kernel of the operator L a — X N in three different linear spaces of formal power 
series. First set 



(4.20) 



considered as a formal power series in y 1 I N (where ^> a is from (|4.5| )). The Airy 
wave function ip a (x,X) (see (p~6|)) is given by the same formula after expanding 



-l/JV 



at x = as in (|4.7| ). The other possibility is to expand y 1 ^ Ar at x = oo: 



(z + ao^r^E 



fc>0 



-i/N 
k 



-i/N-k'-l \N\k 



(4.21) 



Inserting ( 4.21 ) in ( 4.20 ), we obtain another formal series % a (x, A). Denote by ip 

(i) (i) 

ipa , Xa the images of (p a , ip a , \a under the transformations 

yN^^i/N x ^ £ j X) X l/N 



(?) 
a > 



y 
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respectively (e = e 2m / N ). Then ipi. and xi are obtained by expanding (p& and in 
the corresponding spaces of formal series Ker(L a — X N ) d has bases 

{d k x^}, {%<p$}, {%X$h < k < d - 1, < j < N - 1. 

Our observation is that if KerP has a basis 

f l (x) = Y,al j d k z ^\x,z)\ z=x , 

then the same formula gives a basis of KerP when ifi's are substituted by 99's or x's 
and vice versa. Indeed, this follows from ( |1.13j ) and the fact that P has rational 
coefficients. We complete the proof of Proposition |4.9| noting that while P depends 
rationally on x, Xa are formal series in x 1 ^ and the same argument as in the Bessel 
case gives that KerP has a x-basis of the form ( }4.19| ). □ 



5 Explicit formulae and examples 

In this section we have collected several classes of examples. We wanted at least to 
include all previously known examples (unless by ignorance we miss some of them) 
- see |DGj , Wl| , [Z|, p3| , LP |. We hope that we have elucidated and unified them. 



For monomial transformations we derive formulae expressing the operators L 
and A, solving the bispectral problem, only in terms of the matrix A and the vector 



7 (see Proposition 5.1 below). This explicit expression for A (though possibly of 



high order) to the best of our knowledge is new even for N = 2 (see [ DG [). In 
other examples we illustrate the properties of the operator of minimal order from 
a bispectral algebra: when does its order coincide with the rank of the algebra and 
when this operator is a Darboux transformation of a power of a Bessel operator. 
We also point out that the classical Bessel potentials u{x) = cx~ 2 [DG] can produce 
new solutions of the bispectral problem for any c. 

We describe in detail the polynomial Darboux transformations from (L a — X N ) 2 
where L a is an arbitrary Airy or Bessel operator of order N. We do not want simply 
to show that our procedure of constructing bispectral operators works but to point 
out that the involutions a and b (bi in Airy case) possess some very interesting 
properties which deserve further study. 



5.1 Monomial Darboux transformations of Bessel planes 

Let (3 G and W S Gr MB {0). We use the notation from flPH |L2|) (with 
V = Vp) and from ( |3.2| , |3.7| ). When the Darboux transformation is monomial 

d 



9{z) 



h(z) 



(5.1) 



for some n, d. We shall consider only the case when there are no logarithms in the 
basis ( [2.18 ) of KerP. The general case can be reduced to this one by taking a limit 
in all formulae (see Example [O below). Now KerP has a basis of the form 



dN 



fk{x)=^a ki x 1 \ Q<k<n-1, 



(5.2) 
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such that 

7i — 7j £ Nli \ if a fci a fcj / 0, i / j, (5.3) 

where 7 = is from (|2.4j ). 

Let j4 be the matrix [a^i)- We shall use multi-index notation for subsets / = 
{io < . . . < i n -%} of {1, ... , dN} and 5i from (^T7|). We also put 7/ = {7i}ie/, 

= { a k,ii)o<k, l<n-l 

and 
Let 

-^min be the subset of {1, . . . ,dN} with n elements such that det j4^ mm ^ and 
yZicr 7i be the minimum of all such sums, and set 

Eq. ( |5.3|) implies that these numbers are divisible by N. Finally, for a subset / of 
{1, . . . , dN} denote by 1° its complement. 

In the following proposition we express everything entering (|L2l|, |L2|, ^ 



only in terms of the matrix A and the vector 7. Therefore for each A and f3 S 
satisfying ( |5,3j ) (with 7 = /3 d ) we give an explicit solution to the bispectral problem 
(cf. 



Proposition 5.1 In the above notation the operators and the polynomials from 
(gig , |3~Tt1 ) ore given by the following formulae: 



(a) #(z) = z n , 

(b) /(*) = * dAr -, 

Q = (^det^A 7 L 7j0 _ n5/0 x^) (j^det^Aj^)" 1 . 

(c) g b (z) = z n Z det A 1 A IZ P', 

P b = ^detA 7 A 7 i 7/ (^r /7V . 

(d) / b (z) = z dN ~ n Y, det A 7 A/^, 

Q b = ^2detA I A I (L^/ N L 7i0 _ nSl0 . 

Proof. Note that (c) and (d) follow from (a) and (b) (see the proof of Theorem 3.2). 
To prove (a) we note that 

Wr(/ (x), . . . , f n _i{x),9 p (x, z)) 



Vw{x,z) - 

z n Wr(fo(x),...J n -i{x)) 
^det^Wr^)^/^^) fr j. 



53 det A J Wr(ar») 
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The sum is taken over all re-element subsets I = {iq < i\ < . . . < i n -i} C 
{0, 1, . . . , dN — 1}, x 7/ = {x 7i }j g / and ^i(x, z) are the Bessel wave functions (|2.6| ), 
Using (^) and the simple fact 

En(n— 1) 
, ie/Ti a— ) (5.5) 

we obtain (a). 

To prove (b) we shall apply the involution a directly on the tau-function tw of 
the plane W. Recall that j| |gW | 

^(M) = e^^ T(t "|f 1]) , (5.6) 



where [2 x ] is the vector yz , z 2 /2, . . .). The action of a is given by | Wl ] 

T a w(tl,t2, ■ ■ ■ ,tk, ■ ■ •) = T W (ti, —t 2 , . . . , {-l) k ~ l t k , . . .). 



We shall need the formulae [BHY3 



Tw(t) = Edet^A, (5 ' 7) 

and 

Tl (x) = — Wr (x 71 ) r 7 (x) (5.8) 

where rj(t) is the tau-function corresponding to the wave function \&/(x,z) and 
r(x) = t(x, 0, 0, . . .). Applying a to both sides of ([fT?]) and using (|5.6|) and (|2.6|) we 
obtain 

/ x T l ^ Al ^lT a ( 1 +dN8 I -nS)(x)^a(-y+dN5 I -n5)(x,z) 

m aW {x,z) = — -j- — . (5.9) 

L, det ^ J A/r a(7+(iA r ( 5 / _ n<5 ) (x) 

We compute 

0(7 + dA^/ - n£) = 0(7) + dNSpi - (dN - n)5, (5.10) 

which is a Darboux transformation of 0(7). 
The eqs. ( |5.8| ) and (|5.5| ) imply 

44 = — • (5-H) 

Tj(x) xP-> ^ ' 

To apply ( 5. lip in ( |5.9| ) we have to compute p/o but for 0(7) instead of 7. It is a 
simple exercise to see that 

Pi°( a (l)) = Pi(l) =Pi- 

Using this we obtain 



* ( ^ -dN + n ^ detAlA ^ PlL (am i0 , ri ( v 
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Now Proposition 1.7 gives (b) because 

{Lp)* = (-l) N L aW) for € C N 

and 



a ((«(7))/o) = 7/o -nS l0 . 



□ 



In the following example we consider the case when there are logarithms in the 
basis ( [Hp of KerP. 



Example 5.2 Let d = 2, (5 = (1,1,1), 7 
basis 



(3 2 = (1,1,1,4,4,4) and KerP has a 



fo(x) = x 4 , 

fi(x) = a\x + 2a 2 x Inx, 



Using that In x 



gk x e 
4 



ciqx + aixlnx + C12X In 2 x. 



= o we approximate the above functions with 



fo(x,e) = x , 

e) = aix 1+e + 2a 2 e~ V +e - x 4 ), 
/ 2 (x, e) = a x 1+2e + aie _1 (x 



2x 4+€ + x 4 ) 



Consider the Darboux transformation W(e) of Vg( e ), where (3(e) = (1, 1 + e, 1 + 2e), 
with a basis of the operator P(e) consisting of the functions /^(x, e). After changing 
the basis this corresponds to a matrix (cf. (fT^)) 





1 











(: 







2a 2 
















e(ea + ai) 



We apply (5.7) for Twt e y To make the limit e — > we note that the numerator 
and the denominator depend polynomially on e and that (in the notation of (|5.7|) ) 
both ^2,3,6} an d ^"{2,4,5} tend to one and the same Bessel tau-function. So after 
canceling e 3 and setting e = we obtain that tw is a linear combination of 3 Bessel 
tau-functions: 



Tw 



9Qi' r (-2,i,i,4,4,7) + ^8a 2 {ao - Qi)t(_ 2 -2,1,4,7,7) + 4o 2 r (-2 ,-2,-2,7,7,7) 
9of + 18a2(ao — 01) + 4a2 



(5.12) 



As in the proof of Proposition |5.1| from this formula one can compute the operators 
P, Q, Pb and Q^. It is clear that they also can be obtained by taking the limit 
e — > directly in the corresponding expressions for W(e). □ 

From here to the end of the subsection we shall restrict ourselves to the case when 
(3 d = 7 has different coordinates. We choose the following basis of KerL^ (cf. [ |MZ(| ) 



<S> (k _ 1)d+J (x) := wx^-W , 1 < k < N, l<j<d, 



(5.13) 
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where 



N 



fik,i ■= i, fly •■= (j>k,j-i'Yl(Pi- Pk-(j 1 - 



i=l 



In this basis the action of Lp is quite simple: 

j X, / $(k-l)d+j-l, for 2 < j < d 

L ^-iy+j = \ Q ; f ri = i. 

Let a basis of KerP be 

dN 

fk{x) = ^2a ki $i(x), k = 0, . . . ,n - 1. 



(5.14) 



(5.15) 



i=i 



Example 5.3 Let n 

the form: 



<i, /3j — Pj G A^Z for all j and the matrix A = (a^) has 



(i) 
i 

(i) 



(i) 
o 

(i) 
i 



AN) 
AN) 

l l 

(N) 



t 



(TV) 


(N) 
1 



AN) 



n-1 



t W 
L n-2 



,(1) 



t 



(AT) AN) 



n-1 



(5.16) 



n-2 



The type of the matrix is tantamount to the identities Lpfo = 0, Lpf^i = f^, 
k = 1, . . . ,n — 1. Then KerP is invariant under the action of Lp and by Proposition 
|l~5| the operator L = PLpP -1 is differential of order iV and solves the bispectral 
problem. For a generic f3 6 the spectral algebra has rank A" (i.e. it is C[L]). 
This family can be considered as the most direct generalization of the "even case" 
of J. J. Duistermaat and F. A. Griinbaum [pG|| (see also |[M"Z||). When N = 2 our 



example coincides with it but for N > 2 here we present a completely new class of 
bispectral operators. □ 

In connection with the above example we prove the following proposition. 

Proposition 5.4 Let W G GV#(/3) (/3 € C N -generic) be such that Aw contains an 
operator of order N. Then W is a monomial Darboux transformation of Vp, i.e. 

w eGr MB (p)nGA N l 



(5.17) 



Proof. Proposition |1.5| implies that W G Gtb(/3) belongs to GA N ^ iff 

Lp(KerP) C KerP. 



If we suppose that W GrMB(P) then KerP would contain some elements of the 
form ( f2.20p . The action of Lp on them is easily computed: 



LpD k x ^p(x,e l X) 



X N (D x + N) k ^ p {x,^X). 



37 



Thus the linear space span ^D^Sp(x, £l ^)} 0<k<m can be identified with the space 
of polynomials in D of degree < m, with the action of Lp corresponding to P(D) t— > 
X N P(D + N). It is clear that all the L^-invariant subspaces are of the form 



span|L>^ /3 (x,e i A)| 



o<fc<fc 



for some &o- The corresponding polynomial Darboux transformation is trivial in 
the sense that it leads again to the same plane Vp (the operator P = (Lp — X) k ° 
commutes with Lp). Therefore W G Gtmb(I3)- n 

In the same manner as in Example |5.3| , one can build for arbitrary k rank N 
bispectral algebras with the lowest order of the operators equal to kN . 

It is clear that when the matrix A is not of the form ( 5.16| ) (or a direct sum 



of such matrices) then KerP (given by (|5,15|) ) is not invariant under the action of 
Lp. Proposition implies that in this case the spectral algebra does not contain 
operators of order N. The following example is one of the simplest of this type. 



Example 5.5 Let N = 2, (3 = (fa, fa), fa + fa = 1, d = n = 2. We take KerP 
with a basis ( |5.15 ) where 

" 1 a 




A 





1 b 



for some a, b G C, i.e. 



Then 



f (x) = + a<$> 2 {x) = x^ + 2{pi _ p2 + 2) ^ 1+2 , 

h{x) = d> 3 (x) + b^(x) = x& + — b ^ 9 ^ fe+2 - 



Lpfo{x) = axf 3 \ Lph(x) = bxh 



and KerP is not invariant under Lp when ab ^ 0. The spectral algebra Aw '■ 
PL 2 pC[Lp\P~ l consists of operators of orders 4, 6, 8, 10, ... 

This example is also interesting for the fact that it does not require fa — fa € 22 
The generalization for arbitrary N is obvious. □ 

Another example illustrating Proposition |1.5| is the following one. 

Example 5.6 Let N = 2, (3 = (Pi, fa) G C 2 , fa+ fa = 1, fa - fa e 2Z, d = > 
n = 2. We take KerP with a basis (|5.15| ) where 



A 



A 
10 



Xa + b Xb 
a b 



for some a,b,X G C. 

Then it is easy to see that KerP is invariant under the operator Lp + XLp but 
it is not invariant under any polynomial of Lp of degree < 2. On the other hand 
KerP C KerL^ obviously implies L^ +fc KerP C KerP for k > 0. Therefore the 
spectral algebra Aw is the linear span of the operators 

P (Lp + XLp) P" 1 , PLp +k p-\ k>0. 
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This example is interesting for the fact that (for A / 0) the operator of minimal 
order in the spectral algebra is not a Darboux transformation of a power of Lg, 
although the Darboux transformation is monomial. □ 

In the last example of this subsection we show that for d = n = 1 our results agree 
with those of Q . 

Example 5.7 Let d = n = 1, KerP = C/ , 

N f( \ 

and Pi — pj E iVZ if cuaj 7^ 0. Then 

L = PL B P- 1 =PQ=(J2 1 (E ( d * ~ f ) ) X 

A = PbQb = (* - f ) x 

where Pi = Pi - P min , Pmm = mln^, 

iV 

=n^-^)' D x =xd x . 
i=i 

We have 

6(x) = / (J2 a i xPl ) * > deg 6 = iV + 2(/3 max - /3 min ) 
where /3 max = max/?j. When /o(x) = tx^ 1 + x^ 2 , — Pi = Not, a G Z> 

e(x) = x JV (t + x Afa ) 2 

and we obtain the operator A from [g]. □ 

5.2 Polynomial Darboux transformations 

In this subsection we shall consider the simplest case of polynomial Darboux trans- 



formation of an operator of order N, namely when the polynomial h(z) from (1.27) 
is equal to (z — X N ) 2 for some A 6 C\0. Using the kernels of the operators P, Q*, 
P h and Ql from flL2| [L2|, |J, U), we describe the action of the involutions a 



and b (b\ in the Airy case). The Propositions 5.1C, 5.12 below raise some interesting 
questions and conjectures. 

The Bessel and Airy cases are very similar. We shall consider first the Airy one 
since it is simpler. 
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Let Gr A (a), a G C^ 1 . Set 

h(z) = (z-X N ) 2 , g(z)=f(z) = z N -X N . 
Then Kerh(L a ) has a basis of the form 

I J 0<j<AT-l, fe=0,l 

and KerP has a basis 

fj(x) = t> a (x, e j X) + ad x V a (x, e j X), 0<j<N-l 



(5.18) 



(5.19) 



(5.20) 



for some a G C. 

We shall start with the case N = 2. The following example is due to [G3, LP 
We shall obtain it as the simplest special case of Theorem fO| 



Example 5.8 Let N = 2 and a = (a ) G C 1 . For fixed a ,a G C \ we take the 
basis ( 15301) of KerP: 



fk(x) =ipk(x) +ad x ip k (x), k = 0,l 

where ipk{ x ) = ^a(x, (— l) fc A). Using that 

^/fc = a(a 2; + A 2 )^ fc + d^fc, 

= (aao + + X 2 )ip k + a(a x + X 2 )d x i/j k 

we compute P from 

Wr(/ , /!,¥>) 



(5.21) 



PV9 



Wr(/ ,/i 



1 


a 




a(aoa; + A 2 ) 1 




a«o + c^o^ + A 2 0(00^ + A 2 ) d 2 ip 




1 a 






a(aox + A 2 ) 1 





The result is 
P = dl + 



a 2 ao 



-d x + 



a 2 (aox + A 2 ) 2 — (aox + A 2 ) - aa 



1 — a 2 (ao^ + A 2 ) 1 — a 2 (aox + A 2 ) 

This expression coincides with that given in [G3] if we set 

2 2 8 

an = = — , a = — , A = U. 

u 2 + 3t s 2y 

We compute the operators P, Q and Q* as follows. If we write 

P = dl+pi{x)d x +p (x) 
Q = d 2 x +qi{x)d x + q (x) 

Q* = ®l + qi(x)d x + q (x) 



(5.22) 
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then the identity QP = h(L a ) imply 

qi + Pi = 0, 2p[ + qipi +po + 

and 



-2(a x + A 2 ) 



qi 



-qi, qo = -qi + qo- 



Our observation is that because P*Q* = h(L a r a \) and VP a W = f Q*^a(a)i the 
operator Q* has a basis of the form ( |5.2lD with some b G C instead of a and a(a) 
instead of a. Comparing the above expressions for Q* with ( 5.22 ) we obtain that 
b = —a. By Theorem 4.5 the operator P b also has a basis ( jUg ) with some c instead 
of a and \x instead of A. On the other hand we can compute it directly using eqs. 
( [4.15 , 4.16). Then g\>{z) = 1 — a 2 (z 2 + A 2 ) which on the other hand is up to a 
constant z 2 — fi 2 . This gives 

, 1 - a 2 A 2 



(5.23) 



The other coefficients give a surprising result: c = a. In conclusion, if we denote the 
operator P from ( |5.22| ) with P(a, A) then 

P = P(a,X), Q = P*(-a,X), P b = P(a,/x), Q b = P*(-a,/x) (5.24) 

where p, and A are connected by ( 5.23|) . □ 

The next example is completely analogous to the above one but to the best of our 
knowledge it is new. 



Example 5.9 For N = 3 the Airy operator is 

L a = dl + a 2 d x - 



0£qX, 



a = (aoi 0:2) G C 2 , ao ^ 0. We take P with a basis ( 5.20[ ) (N = 3). Then using the 
eq. (14.81) we compute 



P 



a 3 ao 



a 3 (aox + A 3 ) + (1 + a 2 cti 



dl 



+ 



X 



a 3 a2(ttox + A 3 ) + (1 + a 2 «2)«2 + a 2 ao 

a 3 (aox + A 3 ) + (1 + a 2 a.2) 
a 3 (aox + A 3 ) 2 + a«o(l + a 2 «2)(l + a 2 a 2 )(ao^ + A 3 ) 



a 3 (a>ox + A 3 ) + (1 + a 2 «2) 
A direct computation using Proposition |1.7| , Theorem and QP = h(L a ) leads to 
P = P(a,X), Q = -P*(-a,-X), P b = P(a,/x), Q b = -P*(-a, -fj) (5.25) 
with n given by 

p + X = 5-26 

□ 



The above examples can be generalized for arbitrary N as follows. 
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Proposition 5.10 Denote by P = P(a, A) the operator P with a basis ( 5.20| ), Th 
in the above notation we have 



en 



Q = (-lf P *(-a,-X), P h = P(a,ti), Q h = (-l) N P*(-a,-fi) (5.27) 
with X and \x connected by 

X N + ijl n = P a /(—l/a) (5.28) 
where P a > is the polynomial from ( [4.1[ ). The spectral algebras 

Aw = P{L a -\ N ) 2 C[L Q \p- 1 (5.29) 
A blW = Pb {L a - ^) 2 C[L Q ]P b - 1 (5.30) 

consist of operators of orders 2N, 3N, 4N, . . . 

Proof. Because 



i-ifp* (-i) 7V g* = (L a(a) - (-xyr, * aW (x,z) 



N\2 



irQ*y a(a) (x,z) 



-A) 



iV 



QbPb = (L c 



P* 



P h ^ a {x,z) 
z N — n N 



we see that Q = (—l) N P*(b,—X) and Pb = P{c,n) for some b, c, /x. Using the eq. 
L a (x,d x )^ a (x,s j X) = X N W a (x,e j A) we compute p N (x) = Wr(/ , /i, . . . , /at-i) as 
in the proof of Lemma |2.9| . We obtain 

PN (x) = -(- a ) N (a x + X N - P^i-l/a)). 

Eq. $~17\ ) leads to floT^ ) because g h (z) = const • (z N - /j, n ). Applying ( f>?2$ ) for P h 
instead of P, we obtain 

P a >(-l/a) = P a ,(-l/c). 

Note that the map a i— ► c is an automorphism of CP 1 since it is an involution. The 
only solution of the above equation with this property is c = a. 

To compute (— 1) Q*, we note that its second coefficient is equal to that of P 



which is equal to —p' n (x)/pn(x). This, Proposition [L7| and Lemma 47 imply 

ao a (a0o 
a x + X N -P a ,(-l/a) = a(a) x + (-X) N - P a(a y{-l/b) 

which leads to a polynomial equation for b in terms of a and a. Because a i— > 6 is an 
automorphism of CP 1 we obtain that 5 = —a. 

The eqs. ( |5.29 , 5.30| ) follow from Proposition 5A. □ 

We shall find the analog of Proposition 5.1C in the Bessel case. We use the 
notation from the beginning of the subsection with (3 € instead of a and eq. 
flUD modified as follows (cf. (^20|) ) 



(J = o, 



.JV-1, Dfr 



fj(x) = ^ p (x,e j X) +aD x ^ f5 {x,e j X) 
= xd x ). 



(5.31) 



In the next example we shall study the simplest case N = 2. 
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Example 5.11 For iV = 2, j3 = (1 — u, v) the corresponding Bessel operator is 



Lp = x~ 2 (D x - (1 - u))(D x - v )=d 2 x + 



u(l-u) 



x- 



Using ( p., 38 ) we compute the operator P from 

fo fx V 

D x f D x h D 

D 2 J D 2 x h Bl^p 



Pip 



1 



fo fi 
D x fo D x fi 



The answer is the following. If we set 

2 a + 1 + a?v(\ - v) 



j 2 A 2 



then 



P = 2 7 2 , {p2(x 2 )D 2 x + Pl (x 2 )D x +p (x 2 )} 



X 2 p2(x 2 ) 

withp 2 (x 2 ) = x 2 -/j. 2 , pi{x 2 ) = fi 2 -3x 2 &ndp (x 2 ) = -A 2 x 4 + (2A 2 /j 2 + (a + l)(2a- 
l)a~ 2 )x 2 + ((a + l)a -2 - X 2 fi 2 )fi 2 . The operator P b is (cf. 



= -}—{D 2 x p 2 (Lp) + D xPl {Lp) + Po (Lp)\ 

and 5b (-2) = ^ 2 (^ 2 — /J 2 )- A straightforward computation shows that if we set 
P = P(a, A, n) then 

Q = P*(-a/{a + l),\, l i), P h = P(a,n,\)Lp, Q h = LpP*(-a/ (a + 1), /x, A). 

(5.32) 

Therefore we can take 

P h = P(a,»,\), Qb = P*(-o/(o + l),//,A) (5.33) 

i.e. the involution 6 acts simply by exchanging A with /x and vice versa, while the 
involution a acts as a i— > —a/ (a + 1). □ 

The action of the involutions for arbitrary iV is given in the next proposition. 

Proposition 5.12 Denote by P = P(a, A) t/ie operator P with a basis ( 5.31| ) . Then 
we can take P^ and Q\> such that 

Q = (-l) N P*(b,-\), P h = P(a,fx), Q h = (-l) N P*(b,-n) (5.34) 

with A, [a and a, b connected by 

, JV..AT 



A A // x — Pp{—l/a), i + i+AT-l = 



where Pp is the polynomial from ( |1,32| ) . The spectral algeb 

A w = P{Lp-\ N ) 2 C[Lp\p- 1 
A bW = P h {Lp- f i N ) 2 C[Lp]p- 1 
consist of operators of orders 2N, 3N, 4N, . . . 



ras 



(5.35) 

(5.36) 
(5.37) 



43 



Proof. We have g\>{z) = const ■ z (z — fj, ) for some fi G C. Using ( ^ ) we compute 
g h {z) = z 7V det(^/ J (z)) JJ=0 ,...,iv-i = (-a) N z N (z N \ N - Pp(-l/a)) 



which gives the value of [i. We have to prove that P^ given by ( |3.4j ) (which is of 
order 2N) is divisible by Lp from the right. Indeed, it is easy to see that 



p h ( x ,d x )x^ = p{x,d x )x l3i 



A=0 



and the proof of Lemma implies P|a=o = Lp. Thus we can take Pb = P(c, fj) for 
some c G C. Now ( 5.35| ) implies Pp{ — 1/a) = Pp(— 1/c) leading to c = a. 



Finally, as in the Airy case, if Q = (— 1) P*(b, —A) for some b G C then 

°a(/3)( 



P„(-l/a) = (-^^^(-l/ft) = P^l/b + iV - 1) 



showing that a 1 + b 1 + N — 1 = 0. 

The eqs. ( |5.36 , 5.37| ) follow from Proposition 5.4. □ 



In conclusion we want to make some comments. In the case N = 1 the adjoint 



involution a has a simple and beautiful geometric interpretation (see |W1[| ): in terms 
of Krichever's construction it preserves the spectral curve and maps the "sheaf of 



eigenfunctions" into some kind of a dual sheaf. In [Wl] G. Wilson also posed the 
problem of describing the action of the bispectral involution on Gr ad . 

We think that in the general case the study of the action of the involutions a 
and b on the bispectral manifolds of polynomial Darboux transformations of Bessel 
and Airy planes is equally interesting and difficult task. 

The above examples lead us to the conjecture that the involutions a and b (bi 
in the Airy case) possess some universality property. Any polynomial Darboux 
transformation W of a Bessel plane Vn (respectively an Airy plane V a ) is deter min ed 



by the points Ai, . . . , Xn 0) at which the c onditions C are supported (see (2.1)), 
by the matrix A defined by ( fOol) (resp. (|2.19|) ), and of course by the vector (3 (resp. 
a). Then the corresponding matrices for aW and bW (resp. b\W) depend only on 
the matrix A. The point is that they do not depend on the points Ai, . . . 3 Ajy at 
which the conditions C are supported nor on the vector (3 (resp. a). 
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